use "Fig12data.dta" 

gen treat = 1 if total_intensity > .6
replace treat = 2 if total_intensity > .5 & treat == .
replace treat = 3 if total_intensity > .4 & treat == .


. lgraph total_soil_base year if year < 1970, by(treat)

lgraph total_soilbase_area year if year < 1970, by(treat) legend( label (1 "High Crop Intensity") label (2 "Medium Crop Intensity") label (3 "Low Crop Intensity"))

